DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA 93943 



NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




THESIS 



EXPERIMENTAL INVESTIGATION OF TURBULENT 
HEAT TRANSFER IN STRAIGHT AND CURVED 
RECTANGULAR DUCTS 


by 




Joel L. Wilson 




December 1984 




Thesis Advisor: 


M. D. Kelleher 



Approved for public release; distribution unlimited. 

T223280 




Unclassified 



SECURITY CLASSIFICATION of THIS PAGE (Whan Data Entered) 



REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


1. REPORT NUMBER 


2. GOVT ACCESSION NO. 


3. RECIPIENT’S CATALOG NUMBER 


4. TITLE (end Subtitle) 

Experimental Investigation of Turbulent 
Heat Transfer in Straight and Curved 
Rectangular Ducts 


s. TYPE OF REPORT 4 PERIOD COVERED 

Master's Thesis 
December 1984 


6. PERFORMING ORG. REPORT NUMBER 


7. AUTHORfs; 

Joel L. Wilson 

\ 


e. contract or grant number^#; 


9. performing organization name and address 

Naval Postgraduate School 
Monterey, California S3943 


10. PROGRAM ELEMENT. PROJECT, TASK 
AREA 6 WORK UNIT NUMBERS 


II. CONTROLLING OFFICE NAME AND ADORESS 

Naval Postgraduate School 
Monterey, California 93943 


12. REPORT DATE 

December 1984 


13. NUMBER OF PAGES 

81 


14 MONITORING AGENCY name & AOORESSfl/ different from Controlling Office) 


IS. SECURITY CL ASS (of ihlt tapott) 

Unclassified 


1 5a. DECLASSIFICATION DOWNGRADING 
SCH EDULE 


16. Distribution statement (oi i/i/s Rapom; 

Approved for public release; distribution unlimited. 


17. DISTRIBUTION ST ATEmEnT (of the abstract entered In Block 20, If different from Report) 


10 . supplementary notes 


19. KEY WORDS (Continue on reverse side If necessary and Identity by block number) 

Taylor-Gortler Vortices, heat transfer, secondary flow, 
rectangular curved channel, rectangular straight channel, 
Temsheet, Joulean heating, turbulent flow. 


20. ABSTRACT (Continue on reverse side If necessary and Identify by block number ) 

An experimental investigation has been conducted to 
examine the convective heat transfer in straight and curved 
ducts of rectangular cross section. The experimental con- 
figuration was modeled as infinite parallel plates with one 
wall at a constant heat flux and the opposite wall adia- 
batic, Experiments were conducted at steady state for 



Unclassified 



DD F0RM 

\ JAN 73 



1473 EDITION OF I NOV 65 IS OBSOLETE ^ 

S N 0102* LF- OM- 6601 



SECURITY CLASSIFICATION OF THIS PAGE (When Dete Entered) 



Unclassified 



SECURITY CLASSIFICATION OF THIS PACE (TWimi D«<« Enttfd) 



turbulent flow using average Nusselt numbers to compare the 
heat transfer characteristics of the straight and curved 
sections . 

The development of Taylor-Gortler vortices in the 
curved section proved to enhance the heat transfer rate over 
that of the straight section. Improved heat exchanger and 
turbine blade design could both result from a better under- 
standing of the effects of curvature on the rate of heat 
transfer . 



$ N 0102- LF. Ola- 6601 



2 



Unclassified 



SECURITY CLASSIFICATION OF THIS PAGEf»Tl»n />•<• Enl»rmd) 



Approved for public release; distribution unlimited. 



Experimental Investigation of Turbulent Heat Transfer 
in Straight and Curved Rectangular Ducts 



by 



Joel L. Wilson 

Lieutenant, United States Navy 
B.S., United States Naval Academy, 1977 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN MECHANICAL ENGINEERING 

from the 

NAVAL POSTGRADUATE SCHOOL 
December 1984 



John N. Dyei/, Dear/ of Science and Engineering 



3 



ABSTRACT 



An experimental investigation has been conducted to 
examine the convective heat transfer in straight and curved 
ducts of rectangular cross section. The experimental con- 
figuration was modeled as infinite parallel plates with one 
wall at a constant heat flux and the opposite wall adia- 
batic. Experiments were conducted steady state for 
turbulent flow using average Nusselt numbers to compare the 
heat transfer characteristics of the straight and curved 
sections . 

The development of Taylor-Gortler vortices in the curved 
section proved to enhance the heat transfer rate over that 
of the straight section. Improved heat exchanger and turbine 
blade design could both result from a better understanding 
of the effects of curvature on the rate of heat transfer. 
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INTRODUCTION 



I . 

A. TAYLOR-GORTLER VORTICES 

Since the early part of the twentieth century, consider- 
able research has shown that fully developed laminar flow 
along a concave wall does not remain two-dimensional [Refs. 

1, 2, 3] . The flow instead forms a system of spiral vortices, 
of counter rotating pairs, whose axes are aligned in the 
principle flow direction. This phenomenon is the result of 
the variations in the centrifugal forces acting on the 
fluid particles, and is known as Taylor-Gortler vortices. 

Figure 1 illustrates the type of fluid motion just 
described . 

In a channel that is curved in the streamwise direction, 
those fluid particles located in the center of the flow cross- 
section are subjected to higher centrifugal forces than those 
fluid particles traveling along the channel's boundary wall. 

As a result, the fluid in the center of the channel moves out- 
wardly towards the concave boundary. As the process continues, 
the fluid particles near the boundary wall, move in a spanwise 
direction, and finally radially inward replacing the outwardly 
moving particles. These particles then come under the same 
centrifugal forces and the process repeats itself. The re- 
sulting cyclic motion causes the formation of the counter 
rotating Taylor-Gortler vortices, considered to ba primar y 
a laminar flow phenomenon. 
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Figure 1. Schematic of Taylor -Gor tier Vortices in a Curved 
Channel . 
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It has been observed [Ref. 4] that the heat transfer rate 
from flow along a concave curved wall is greater than that 
for flow along a straight wall, a phenomenon attributed to 
the additional mixing provided by the secondary motion of 
the Taylor-Gortler vortices. It is thought that the cross- 
hatching observed on reentry vehicles can be explained in 
part by the presence of streamwise vortices similar to the 
Taylor-Gortler vortice phenomenon. 

There are many possible applications that could result 
from a more thorough understanding of Taylor-Gortler vortices 
and their effect on heat transfer and fluid flow character- 
istics. Two such applications could lead to improved heat 
exchanger designs and improved turbine blade cooling [Refs. 

5, 6, 7] . 

B. HISTORY 

The instability of an in' iscid fluid flowing past a 
curved boundary, was first considered by Lord Rayleigh in 
1916 [Ref. 8]. By assuming that the fluid was non-viscous, 
he determined that for the motion of an inviscid fluid to 
remain stable, its circulation must increase with increasing 
radius. G. I. Taylor, in 1923, [Refs. 1, 9], continued this 
study with an extensive analytical and experimental study 
of viscous fluids. His investigations focused on the flow 
between two cylinders, in which the inner cylinder rotated 
while the outer cylinder remained stationary. Taylor 
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ascertained that such flows become unstable when the value 
of a dimensionless parameter exceeded a critical value of 
41.3. The parameter, known as the Taylor number is defined 
as : 

Ta = Re 

where 'd' is the width of the gap, assumed small when compared 
to ' Ri ' , the radius of the inner cylinder, and 'Re' is the 
Reynolds number based on the peripheral velocity of the inner 
cylinder. Taylor determined that for those cases in which 
the value of the Taylor number exceeded the critical value, 
a secondary motion developed and the Taylor vortices formed. 
Figure 2 illustrates this fluid motion. 

Instability of a similar nature is also observed when a 
viscous fluid flows in a curved channel due to a pressure 
gradient acting along the channel wall. This problem was 
first considered analytically by W. R. Dean [Ref. 10] in 1928, 
for a channel formed by two concentric cylinders, where the 
radius of the inner cylinder was large in comparison to the 
spacing between the inner and outer cylinder walls. Dean 
concluded that there would be an initiation of flow instabil- 
ity, and the formation of vortices, when a dimensionless 
parameter, the Dean number, exceeded a value of 36. The 
Dean number is defined as: 

„ ,~d~ 

De = Re / 
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Figure 2. Schematic of Taylor Vortices between Cylinders. 



where ' d' represents the channel half-width, ' Ri ' is the 
inner cylinder radius, and 'Re' is the Reynolds number based 
on the mean velocity of the undisturbed flow and the channel 
half width 'd'. The analytical work of Dean was later veri- 
fied by W. H. Reid [Ref. 11], using an approximate solution. 

In 1940, H. Gortler [Ref. 2] studied the stability of 
laminar boundary layer profiles on curved walls under the 
influence of small disturbances. He found that these dis- 
turbances were similar to the vortices studied by G. I. 
Taylor. Using approximate numerical calculations, Gortler 
concluded that the disturbances or vortices were produced 
only on the concave boundary walls and that the overall flow 
profile appeared to remain laminar in nature. These results 
were verified, with an exact solution, by G. Hammerlin, and 
reported by H. Schlichting [Ref. 12] in 1955. A. M. 0. Smith 
completed an even more extensive numerical analysis that fur- 
ther substantiated these findings [Ref. 3], The results of 
these numerical solutions have recently been demonstrated 
with the use of hot wire anemometry, laser doppler systems, 
and flow visualization techniques, [Refs. 13, 14]; and in 
1976, Y. Aihara [Ref. 15] conducted a non-linear analysis 
of Gortler vortices. 

As interest began to develop concerning the effects of 
the secondary flows associated with the Taylor-Gortler 
vortices, many studies were published concerning the in- 
fluences of these vortices on heat transfer in the laminar 
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and turbulent flow regimes. In 1955, F. Kreith [Ref. 16], 
studied the influence of heat transfer with respect to the 
curvature of the boundary wall for fully turbulent flows. 

He concluded that the heat transfer from the heated concave 
boundary wall was considerably higher than that transferred 
from the convex boundary wall of the same curvature and 
under similar turbulent conditions. 

In 1965, L- . Persen [Ref. 17], considered the special cases 
of very igh and very low Prandtl number fluids and related 
the inc se in heat transfer from a curved wall to the pre- 
sence of the Taylor-Gortler vortices. 

There has been only a limited amount of published litera- 
ture dealing with the flow and heat transfer in curved channels 
of rectangular cross-section and large aspect ratios, where 
the aspect ratio is equal to the spanwise distance of the 
channel divided by the channel height. Much of what has been 
published involves the development of numerical approximations 
and solutions for heat and mass transfer in curved ducts of 
various geometries. K. Cheng and M. Akiyama [Ref. 18] devel- 
oped a numerical solution for forced convection heat transfer 
with laminar flows in curved channels of rectangular cross- 
section, but only for small aspect ratios. 

In 1976, A. A. Shibani and M. N. Ozisik [Ref. 19], using 
matched asymptotic expansion techniques for a wide range of 
Prandtl nv )ers, solved the heat transfer problem between paral- 
lel plates vith turbulent flow, for the case of uniform wall 



temperature . 
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Other experimental and analytical studies that are worth 
citing with regard to this present study follow. Y. Mori 
[Ref. 20], obtained results for hydrodynamically fully devel- 
oped flows with constant wall heat flux, in curved channels 
of square cross-section. W. M. Kays and E. Y. Leung [Ref. 21] 
reported solutions for turbulent flow heat transfer in a con- 
centric circular tube annulus with a fully developed velocity 
profile and constant heat rate per unit length, for a fluid 
of Prandtl number 0.7. Results for large aspect ratio 
channels with rectangular cross-sections for laminar flows 
were reported by M. Durao [Ref. 22] and J. Ballard [Ref. 23]. 
R. Holihan, Jr. [Ref. 24] reported results for laminar and 
transition flows and S. F. Daughety [Ref. 25] reported 
results for turbulent flows. P. F. Brinich and R. W. Graham 
[Ref. 26] reported results for turbulent flows in a rectan- 
gular curved channel with an aspect ratio of 6, for the 
inner wall heated, the outer wall heated, and both walls 
heated . 
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II. INTENT C. THE STUDY 



The purpose of this investigation was to examine the 
effect of streamwise curvature on the heat transfer rate in 
a curved rectangular duct of large aspect ratio, again defin- 
ing aspect ratio as the spanwise distance divided by the 
channel height. The flow velocities examined were in the 
transitior and turbulent regime. The results of this study 
were compared to the heat transfer rates for the same range 
of velocities, in a straight duct of identical aspect ratio. 

Although Taylor-Gortler vortices are considered to be a 
laminar flow phenomenon, enhanced heat transfer in curved 
channels for transition and turbulent flows has been observed. 
Even though it has not been proven that the vortices continue 
to propagate at the higher flow velocities, this enhancement 
in the transfer of heat has been attributed to the secondary 
flow velocity components of the Taylor-Gortler vortices. It 
is believed that these secondary components transport the 
heated fluid from the outer concave wall, inward toward the 
opposite wall of the channel, displacing the cooler fluid 
particles and causing them to move toward the heated concave 
wall. It was expected that similar results would be observed 
in this study. 

This investigation was conducted using a single channel 
w_th a rectangular cross-section and constant aspect ratio. 

The channel incorporated both a straight test section and a 
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curved test section. The results obtained in each test section 
at approximately the same flow rates were compared in an effort 
to determine the effects of the Taylor-Gortler vortices on the 
transfer of heat. Also, the results of both the straight and 
curved sections were compared to the results of Ballard [Ref. 
23], Holihan [Ref. 24], and Daughety [Ref. 25]. 

The straight section results of this study were compared 
to the Dittus-Boelter equation for flows in circular tubes 
[Ref. 26], the results of Kays and Leung for turbulent and 
transition flow in annular passages [Ref. 21] , and the 
analytical solutions of Petukhov and Popov for transition 
flows in circular tubes [Ref. 27] . 

The curved section results of this study were compared 
with the results of Brinich and Graham [Ref. 28] for turbulent 
flow in a rectangular curved channel, and also to the solutions 
of Petukhov and Popov, mentioned above. 
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III. EXPERIMENTAL WORK 



A. DESCRIPTION OF THE APPARATUS 

To meet the objectives of this investigation, a channel 
of rectangular cross-section was used. Details of its con- 
struction are described in references 22 and 23. Plexiglas, 
0.635 centimeters thick, formed the inner and outer walls of 
the channel. These walls were separated by 0.635 centimeter 
spacers that also served as the sides of the channel. 

The channel, shown in Figure 3, had a straight section, 
122.0 centimeters in length, followed by a curved section of 
180 degrees of arc, followed by another straight section, 91.4 
centimeters in length. The inner curved wall of the channel 
had a radius of curvature of 30.5 centimeters. The rectan- 
gular channel was 0.635 centimeters high and 25.4 centimeters 
wide resulting in an aspect ratio of 40. The hydraulic dia- 
meter was 1.229 centimeters. The cross-sectional area of 
the channel was 16.13 square centimeters. The wetted peri- 
meter was 52.07 centimeters. Figure 4 shows a cross-sectional 
view representing both the straight and curved sections of 
the channel . 

An entrance bell constructed of plexiglas was connected 
to the straight section of the channel. It was designed and 
manufactured according to ASME nozzle standards, with an 
elliptical curved base on a major axis equal to ten inches 
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Figure 3. Illustration of Equipment and Apparatus. 
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and a minor axis of one inch. Cheese-cloth was attached to 
the entrance nozzle to prevent foreign matter from entering 
the flow channel. An aluminum exhaust nozzle was attached 
to the exit of the channel and directed the flow from the 
channel into two inch diameter pvc piping. The two inch pvc 
piping contained a measuring section which could be fitted 
with standard ASME concentric orifices. In this study two 
orifice plates were used. An orifice with a diameter of 
0.5325 inches was used for hydraulic Reynold's numbers 
below 5000, and an orifice with a diameter of 1.065 inches 
was used for hydraulic Reynold's numbers above 5000. Pressure 
taps on either side of the orifice, at 1 and 1/2 pipe diameters, 
and connected to manometers were used to calculate the mass 
flow rate of the working fluid, which was air at room 
temperature . 

The flow of air was drawn through the channel by an 
electrically driven Spencer Turbo Compressor, rated 30 horse- 
power at 3500 rpm, and 550 cubic feet per minute at 70 degrees 
fahrenheit and one atmosphere. 

The channel contained two test sections from which experi- 
mental heat transfer data could be obtained. The straight 
test section was located 120 centimeters downstream of the 
entrance bell to ensure hydrodynamically fully developed 
flow. The straight test section was 29.2 centimeters in 
length with a heated test section area of 741.7 square centi- 
meters. The curved test section was located in the lower 
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half of the curved portion of the channel. It was 28.3 
centimeters in length, subtending an arc of 53.1 degrees, 
and having a heated test section area of 718.4 square 
centimeters . 

In each of the heated test sections, Temsheet, a carbon 
impregnated porous paper with a uniform electrical resist- 
ivity, was glued to the outer wall. Joulean heating was 
used to heat the flow of air through the channel. Since the 
electrical resistance of the Temsheet is not constant, but 
varies slowly with temperature, a precision resistor with an 
electrical resistance of 2.0262 ohms was connected in series 
with the Temsheet to allow the calculation of the instantan- 
eous power being supplied. Detailed schematics of each test 
section are shown in Figures 5 and 6. 

The variables that were measured and used in this invest- 
igation were: 

(1) the temperature of the air entering the channel (T^ ) 

(2) the temperature of the air leaving the channel (T ) 

(3) the temperature of the heated outer wall for each 

test section (T ) 
wo 

(4) the temperature of the unheated inner wall for each 
test section (T . ) 

W 1 

(5) the temperature between the three layers of insulation 

at each test section (T. ) 

ins 

(6) the voltage across the precision resistor (V ) 

(7) the voltage across the Temsheet or heater (V^) 
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Figure 5. Straight Test Section, Detailed Schematic. 
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Figure 6. Curved Test Section, Detailed Schematic, 
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(8) the temperature of the air at the orifice ( T or jf) 

(9) the pressure upstream of the orifice (P ) 

(10) the difference in the upstream and the downstream 
pressure across the orifice (AP) 

(11) the atmospheric pressure (P a ^_ m ) 

A total of sixty-one Copper-Constantan , 30 gauge thermo- 
couples were located throughout the test apparatus to record 
the desired temperatures. Readings from four thermocouples 
spanning the width of the channel were averaged to give the 
temperature at the entrance and exit to the channel. Five 
thermocouples were connected in parallel and inserted between 
each of the layers of insulation at each of the test sections. 
In addition, individual thermocouples were positioned in each 
of the test sections, eight in direct contact with the Temsheet 
and eight in the plexiglas of the unheated inner wall. The 
beads of the thermocouples in contact with the Temsheet were 
electrically insulated with ENMAR Heat Resisting Glyceryl 
Phthalate. Small diameter holes were drilled in the plexi- 
glas to allow the installation of the thermocouples. One 
thermocouple was inserted in the two inch pvc pipe downstream 
of the orifice in accordance with the ASME Power Test Code 
to record the temperature of the air flowing through the 
orifice . 

The thermal insulation consisted of 1/2 inch layers of 
Armstrong Armaflex 22 Sheet Insulation, a flexible foamed 
plastic material. By positioning thermocouples between 
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these insulation 1 yers, the eat lost to the environment 
could be computed t each tes_ section. The heat loss 
through insulation surrounding the rest of the channel due 
to the heated air flowing through the channel was small in 
comparison. See Figure 7 for a detailed sketch of the 
thermocouple placement and attachment. The insulation was 
held in place by the use of velcro strips attached to the 
frame. The entire channel and all connections between the 
channel, the pvc piping, and the orifice were sealed with 
General Electrical Silicone Rubber Sealant Caulk to ensure 
there was no leakage of air into the channel or the piping 
which could affect the recorded temperatures or mass flow 
rate . 

The data acquisition system used for this experiment 
was a Hewlett Packard 3054A Automatic Data Acquisition/ 
Control System consisting of a 3456A Digital Voltmeter and 
a 3497A Data Acquisition/Control Unit. Also used in con- 
junction with this data acquisition system, were a Hewlett 
Packard 9326 Computer terminal and 2671G Printer. Pressure 
measurements were taken from Meriam vertical manometers, one 
whose fluid was water with a 0 to 60 inch range, another 
using water with a 0 to 30 inch range, and a third whose 
fluid was mercury, calibrated to read inches of water with 
a range of 0 to ^15 inches. A photograph of the channel 
and associated st equipment are shown in Figures 8, 9, 
and 10. 
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Figure 7. Thermocouple Placement in Test Sections. 
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Figure 8. Photograph of Test Channel. 
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Figure 10.' Photograph of Data Acquisition System. 








Sheets of cholesteric liquid crystal slurries with a 
temperature range of 40-43°C were used to provide qualita- 
tive measurements of the surface temperature distribution 
on the Temsheet. The liquid crystals exhibit dramatic color 
changes over a small temperature differential and allow a 
reasonable comparison to be made between the heat transfer 
processes on each test section by observation of the color 
pattern displayed for varying temperature gradients and flow 
velocities. Specifically, the occurrence of multi-colored 
stripes aligned in the flow direction on the curved wall 
would indicate regions of high and low temperatures and con- 
firm the presence of Tay lor-Gortler vortices. 

The addition of this liquid crystal sheet necessitated 
a calculation of its overall heat transfer coefficient. To 
determine the true surface temperature of the heated wall 
the temperature drop across the sheet must be known. The 
sheet was composed of 4.5 mils of a liquid crystal slurry 
with approximately the same thermal conductivity as body 
fat, .5 mils of acrylic backing, and a 2 mil mylar cover; 
resulting in an overall heat transfer coefficient of 
522 w/°C. Insulation opposite the heated sections was 
fitted with removable panels to allow viewing and photo- 
graphing of the liquid crystals. 

B. EXPERIMENTAL PROCEDURES 

The experimental procedures followed were the same for 
the straight and curved test sections. The results obtained 
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from the straight test section served as a baseline for the 
comparison of the curved section results. Experiments were 
conducted by varying the mass flow rate of air corresponding 
to Reynolds numbers from approximately 800 to 25000, where 
the Reynolds numbers are based on hydraulic diameter. For 
each mass flow rate, data was taken automatically by the data 
acquisition system and the data was immediately reduced and 
printed. Preliminary runs were performed to determine the 
time required for the test rig to come to steady state. 

For this experiment three hours were allowed for steady 
state to be reached. After three hours, data was taken at 
ten minute intervals to ensure that a steady state condition 
had been achieved. The criteria for steady state was based 
on three variables: 

(1) the difference between inlet and outlet temperature 

(2) the heated boundary wall temperature 

(3) the unaccounted heat loss in the system 

When the first two of these variables varied by less than 
two percent over a ten minute interval and the third began 
to fluctuate, it was considered that a steady state condition 
had been achieved. 

It was also necessary to determine the approximate power 
setting required to bring the heated wall of each test section 
to approximately 50 degrees Centigrade. It had been deter- 
mined by Holihan [Ref. 24] that this temperature was suffi- 
cient to ensure a 20 degree Centigrade difference in the 
heated wall and the unheated wall. 
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For each run, the atmospheric pressure, the pressure 
difference across the orifice, and the pressure upstream of 
the orifice were entered into the computer for the mass flow 
rate calculations. All other data was acquired by the data 
acquisition system directly from the experimental apparatus. 

To compute the instantaneous power supplied the following 
relationship was used: 



P r pr 



The precision resistor voltage and the heater voltage were 
read by the data acquisition system and the resistance of the 
precision resistor w as a known constant. A diagram of this 
circuit is shown in Figure 11. 

Upon completion of selected data runs photographs were 
taken of the liquid crystals on the heated section to provide 
a visual record of the airflow patterns. 

For Reynolds numbers ranging from approximately 2400 to 
5600, it was necessary to make two sets of data runs for the 
straight test section. The first, with conditions the same 
as for other flow rates, which exhibited a pronounced devia- 
tion from the near linear results; and a second, where the 
flow just prior to the straight test section was tripped 
with sandpaper, which reestablished the near linear results. 
This reproducible anamoly indicated the presence of a tran- 
sition region between laminar and turbulent flow. 
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Power Supply 



Figure 11. Circuit Used to Measure Power Supplied to 
the Temsheet. 
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IV. PRESENTATION OF DATA 



A. ANALYSIS 

The constant heat flux surface for the heated wall of each 
test section was provided by the uniform electrical resis- 
tivity of the Temsheet. The insulated unheated wall was con- 
sidered adiabatic since the heat losses through that wall were 
negligible. The channel's large aspect ratio of 40 allowed 
the modeling of the channel as parallel plates. The straight 
portion of the channel upstream of the straight test section 
was of sufficient length to ensure that the flow was hydro- 
dynamically developed for the flow velocities of this study. 
The straight portion of the channel downstream of the curved 
test section was also of sufficient length to ensure that the 
flow exiting the channel was thoroughly mixed and that the 
temperature, T , was an average bulk temperature. 

Based on these assumptions, the experimental configura- 
tion was modeled as forced convection between parallel plates, 
for hydrodynamically developed and thermally developing flows, 
subjected to a constant heat flux. The boundary conditions 
are one wall at a constant heat flux and the other wall 
insulated . 

To analyze this problem, several quantities were defined 
as follows; 

The heat convected to the air was calculated using the 



equation : 
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■air 



= m C 



pa 



(T , - T . ) 

.r out in 



where 1 C . ' was the specific heat of the air at constant 

pci 1 IT 

pressure and 'm' was the mass flow rate of the air. Mass 
flow rate was calculated using standard ASME flow measure- 
ment techniques [Ref. 33]. 

The average heat transfer coefficient between the heated 
wall and the flow of air in the channel was def .ed by the 
equation : 



h = 



^air 



a pl at 



where 'Q ■ 'is defined above, 'A^ T ' was the area of the 
air PL 

Temsheet in the test section, and 'AT' was the difference 
between the arithmetic average heated wall temperature (T ) 
and the average bulk air temperature (T^p. The average 
bulk temperature was defined as the arithmetic mean of the 
entrance and exit temperatures (T^ , T Qut ) • 

The average Nusselt number was then calculated using: 



Nu 



hd 



h D 



hd 



K 



air 



In this equation ' D, is the hydraulic diameter and 'K . ' 

n a ^ i y* 



air 



is the thermal conductivity of air. 

The Reynolds number was calculated for each test run as 
follows : 



m D 



Re 



hd 



hd 



A y . 
c air 



42 



where again 'm' and ' 1 are the mass flow rate and hydrau- 
lic diameter of the channel, ' A ' is the cross-sectional 
area of the channel, and ' u • 'is the dynamic viscosity of 

3 . 1 1C 

the air. 

For the curved section runs the Dean number was defined as: 



where ' Re^' is the Reynolds number based on the hydraulic 
diameter, ' ' is the hydraulic diameter, and ' Ri ' is the 
radius of curvature of the unheated inner convex wall. 

A sketch of the control volume and a set of sample calcu- 
lations for one test run of the curved section are given in 
Appendix B. A sample calculation for the uncertainty analysis 
is given in Appendix A. 

B. RESULTS 

The data obtained from each experimental run was evaluated 
utilizing the expressions described in Part A above. The 
major parameters resulting from this evaluation are shown in 
Tables I through V. Tables I through III contain the straight 
section results and Tables IV and V contain the results of 
the curved section. A plot of the average Nusselt number 
verses Reynolds number is given in Figure 12 for the compari- 
son of the results. Figure 12 shows the flow passing through 
the transition regime of the straight section in both the 



De = Re 



hd 
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Table I 



Summary of Straight Test Section Results 
(Large Orifice) 



Re,, 

hd 


^air 

(W) 


h 

(W/mC) 


AT 

(C) 


Nu hd 


4973 


47.09 


30.80 


20.61 


14.97 


5633 


67.70 


35.65 


25.60 


17 .26 


6326 


83.83 


39.19 


28.83 


18.81 


7096 


75.40 


44.41 


22.88 


21.54 


7989 


98.82 


47.99 


27.75 


23.09 


8935 


82.48 


52.05 


21.36 


25.10 


9969 


101.85 


56.20 


24.43 


27.29 


11520 


121.29 


63.12 


25.90 


30.52 


12772 


120.00 


67.01 


24.13 


32.35 


14357 


122.20 


73.52 


22 . 40 


35.43 


16137 


122.09 


77.78 


21.16 


37 . 50 


17990 


124.29 


83.17 


20.14 


40.55 


20629 


125.30 


89 . 53 


18.86 


43.82 


24507 


127.56 


101.13 


17.00 


48.63 



Reynolds number and Nusslet number are based on hydraulic 
diameter . 
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Table II 



Summary of Straight Test Section Results 
(Small Orifice - Untripped) 



Re hd 


Q • 
air 

(W) 


h 

(W/mC) 


AT 

(C) 


Nu hd 


1524 


27.17 


14.74 


24.84 


7.20 


2162 


32.19 


18.92 


22.93 


9.20 


2686 


34.81 


19.03 


24.65 


9 . 23 


3321 


39.19 


21.67 


24.37 


10.57 


4063 


40.31 


23.27 


23.34 


11 . 00 


4463 


39.28 


23.83 


22.21 


11.60 


4851 


52.09 


26.41 


26.58 


12.77 


5223 


51.15 


30.97 


22.25 


15.18 



Reynolds number and Nusselt number are based on hydraulic 
diameter . 
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Table III 



Summary of Straight Test Section Results 
(Small Orifice - Tripped) 



Re, , 
hd 


Q 

air 

(W) 


h 

(W/mC) 


AT 

(C) 


Nu hd 


1976 


41.74 


17.91 


31.40 


8 . 62 


2354 


34.35 


19.87 


23.30 


9.62 


2793 


40.11 


22.55 


23.97 


10.88 


3225 


35.37 


24.52 


19.44 


11.95 


3626 


50.87 


28.29 


24.23 


13 .72 


4639 


50.47 


33.84 


20.10 


16.44 


5552 


51.99 


39.33 


17.81 


19.18 


7014 


72.53 


45.85 


21.32 


22.27 



Reynolds number and Nusselt number are based on draulic 
diameter . 
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Table IV 



Summary of Curved Test Section Results 
(Large Orifice) 



Re, , 
hd 


De 


Q • 
air 

(W) 


K 

(W/mC) 


AT 

(C) 


Nu hd 


4952 


2257 


79.49 


41.71 


26.11 


20.39 


5650 


2574 


72.36 


44.95 


22.05 


21.97 


6560 


2989 


82.47 


50.05 


22.57 


24.20 


7578 


3453 


81.40 


54.48 


20.47 


26.39 


8444 


3847 


97.95 


59 . 04 


22.73 


28.56 


10000 


4555 


103.57 


64.62 


21.96 


31.11 


10972 


4999 


105.79 


68.81 


21.06 


33.80 


11515 


5256 


108.32 


72.74 


20.40 


34.97 


12067 


5498 


110.84 


72.50 


20.94 


35.41 


13000 


5922 


114 . 52 


79.05 


19.84 


38.19 


13723 


6252 


116.24 


82.86 


19.22 


40.34 


14493 


6602 


118.59 


86.79 


18.72 


42.36 


17238 


7854 


121 .95 


99.76 


16.75 


48.41 


20517 


9348 


124.65 


113.84 


15.00 


55.50 


24208 


11029 


126.31 


126.43 


13.68 


61.20 


Reynolds 
diameter . 


number and 


Nusselt : 


number are based 


on hydraul 


ic 
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Table V 



Summary of Curved Test Section Results 
(Small Orifice) 



Re hd 


De 


Q • 
air 

(W) 


h 

(W/mC) 


AT 

(C) 


Nu hd 


1550 


706 


39.59 


22.75 


23.84 


10.86 


1828 


833 


40.55 


24.86 


22.35 


12.04 


2150 


978 


54.24 


27 . 41 


27.11 


13.10 


2426 


1105 


50.84 


29.66 


23.48 


14.27 


2751 


1254 


48.48 


30.38 


21.86 


14.83 


3372 


1536 


49.48 


33.56 


20.20 


16.00 


3815 


1738 


51.06 


36.39 


19 .22 


17.65 


4312 


1964 


71.12 


39.11 


24.91 


19.05 


5215 


2378 


72.52 


44.02 


22.57 


21.22 


Reynolds 


number and 


Nusselt 


number are based 


on hydraulic 



diameter . 
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1000 Rchd luuou 

Figure 12, Straight vs. Curved Section Results Present Study. 



tripped and untripped condition. All further comparisons 
will be conducted using the data from the tripped flow. 

The results indicate an increase in the rate of heat 
transfer with increasing Reynolds number for both the 
straight and curved test sections. In addition, the heat 
transfer rate was higher in the curved section for each 
Reynolds number investigated. 

A straight line least squares correlation of the data 
for the present study resulted in: 

Nu = 0.02 Re hd 0 ' 79 2300 < Re hd < 7100 

for transition flow in the straight test section; 

Nu = 0.22 Re hd ° " 53 1500 < Re h(J < 5300 

for transition flow in the curved test section; 

Nu = 0.07 Re hd °* 65 8900 < Re hd < 25000 

for turbulent flow in the straight test section; and 

Nu = 0.03 Re hd °' 76 8400 < Re hd < 25000 

for turbulent flow in the curved test section. 

In the short laminar regime of this study the rate of 
heat transfer was about forty-five percent higher in the 
curved section, indicating the presence of Taylor-Gortler 
vortices. As flow progressed through the transition regime 
a loss of enhancement of heat transfer rate, from forty to 
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about twenty percent, was observed as flow over the straight 
section went from laminar to fully turbulent. In the low 
end of the turbulent regime, Reynolds number between 5400 
and 12000, the rate of heat transfer was about thirteen 
percent higher in the curved section. As the Reynolds 
number continued to increase beyond 12000 so did the enhance- 
ment of heat transfer, to approximately twenty-six percent 
at 24000. This compares favorably with earlier studies by 
Kreith [Ref. 16] , where he reported an increase in the rate 
of heat transfer along a concave wall of from twenty-five 

to sixty percent for Reynolds numbers, based on hydraulic 

4 6 

diameter, between 10 and 10 . Holihan reported an increase 
of fifteen percent for laminar flow and thirty percent for 
transition flow [Ref. 24]. Ballard reported an increase of 
eleven percent in the heat transfer rate for laminar flows 
[Ref. 23] . Daughety reported an increase of twenty percent 
in the heat transfer for turbulent flows [Ref. 25] . 

The actual presence of Taylor-Gortler vortices was 
verified by the observation of the liquid crystal isotherm 
distribution in the curved section. At low values of 
Reynolds number and through most of the transition regime, 
liquid crystal bands oriented in the streamwise direction 
were visible (Fig. 13). As the Reynolds number increased 
turbulent forces began to dominate in the curved section, 
and less evidence could be detected for the existence of 
the vortices (Fig. 14) . Over the entire range of Reynolds 
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numbers studied for the straight section, the liquid crystals 
formed essentially straight bands perpendicular to the flow 



direction , 
(Fig. 15). 
due to the 
Figures 16 



indicating approximately two-dimensional flow 
Disruption of flow in the transition region 
insertion of the sandpaper is illustrated in 
and 17. 
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Figure 13. Liquid Crystal Representation of Airflow. 
(Curved Section, Re^ = H00) 




Figure 14. Liquid Crystal Representation of Airflow. 
(Curved Section, Re^ = 9400) 
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Figure 15. Liquid Crystal Representation of Airflow. 
(Straight Section, Re^ d = 12700) 
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Figure 16. Liquid Crystal Representation of Airflow. 

(Straight Section, Re hd = 4000, Untripped) 




Figure 17. Liquid Crystal Representation of Airflow. 

(Straight Section, Re^ = 4000, Tripped) 
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V. DISCUSSION AND CONCLUSIONS 



Holihan [Ref. 24] , in his study, determined that there 
was a negligible difference in the fluid bulk temperature 
and the unheated wall temperature. He additionally demon- 
strated that the radiated heat transfer from the unheated 
wall surface was minimal. These results were verified by 
Daughety [Ref. 25] for the present study and provide the 
basis for the assumption that the heat transfer to the air 
flowing through the channel was solely by convection from 
the heated Temsheet. 

As was mentioned earlier, the high aspect ratio of the 
channel provided the basis for the assumption that the ex- 
perimental apparatus, as configured, could be modeled as 
infinite parallel plates. Holihan's data [Ref. 24] also 
tended to verify this assumption. His experimental data, 
in the laminar flow region, approached the theoretical limit 
for average Nusselt number for parallel plates with one 
wall heated at a constant heat flux and the opposite wall 
adiabatic [Ref. 29]. 

Based on the assumptions mentioned above, comparisons 
were made with other analytical and experimental results 
that were of the same or similar problem. Again the problem 
being, flow between infinite parallel plates with one wall 
at a constant heat flux and the opposite wall adiabatic. 
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A comparison of the experimental results of this study 
and the experimental results of Ballard [Ref. 23], Holihan 
[Ref. 24] and Daughety [Ref. 25] are shown in Figure 18, 
for the straight section, and Figure 19 for the curved 
section. Figure 20 shows a compilation of all the data on 
one plot. Ballard's data was strictly in the laminar 
region; Holihan' s data covered the laminar as well as tran- 
sition regions, and Daughety 's data was strictly in the 
turbulent region. The channel aspect ratio was 40 for each 
of these studies and the experimental procedures were similar 
in each case. Ballard and Holihan compared their data to 
the analytical studies by McCuen, [Ref. 30] , for heat 
transfer between infinite parallel plates with constant 
wall temperatures and heat flux. Daughety compared his 
data with the results of Shibani and Ozisik [Ref. 19] , and 
Kays and Leung [Ref. 21]. Ballard's and Holihan's data 
plotted above analytical solutions while Daughety 's data 
plotted below. Differences were attributed to the variations 
in the geometries used in the studies, side wall effects, 
inability to account for all of the heat transfer processes, 
and the limitation inherent in any experimental work. 

A comparison of analytical and experimental results 
with those of the present study for the straight test sec- 
tion is shown graphically in Figure 21. The results of 
Petuhkov and Popov [Ref. 27] are shown for turbulent heat 
transfer in round tubes. The correlation used was 
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Figure 18. Comparison of Straight Section Present Data with Daughety 
Holihan, and Ballard. 
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Figure 19. Comparison of Curved Section Present Data with Daughety, 
Holihan, and Ballard. 
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200 1000 R^hd 10000 

Figure 20. Comparison of Both Sections Fresent Data with Daughety, 
Holihan, and Ballard. 



Nu 



(f/8) Re Pr 



k 1 (f)+k 2 (Pr) v'f/8 (Pr (2/3) -l) 



where ' f ' the friction factor, has been calculated using the 
Filonenko equation given in Karlekar and Desmond [Ref. 32] as 



and 



f = (1.82 log Re = 1.64) 



k x (f) 
k 2 (Pr) 



1 + 1.4 f 



11.7 



+ 



1.8 

Pr< 1/3) 



The Dittus-Boelter equation [Ref. 26] for heat transfer in 
a straight tube for constant wall temperature is plotted 
from: 

. T r\mn 0.08 

Nu = 0.02 Re 



For a Prandtl number equal to 0.71. Additionally, the 
results of Kays and Leung [Ref. 21] for heat transfer in 
annular passages are shown for r* = 1.0, which equates to 
parallel plates. When extrapolated down these results 
continue to correlate well with present data. 

The curved section data is plotted in Figure 22, and 
is compared to the experimental work of Brinich and Graham 
[Ref. 28]. Also plotted as a reference are the analytical 
solutions of Petukhov and Popov [Ref. 27] . The accuracy of 
the data points used from Brinich and Graham are subject to 
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Figure 21. Comparison of Straight Section Present Data with Di ttus-Boelter 
Petukhov and Popov, Kays and Leung, an Extrapolated Kays and 
Leung, and a Calculated Least Squares Fit. 



BSEmOOUlCED AT GOVERNMENT EXPENSE 




63 



Figure 22. Comparison of Curved Section Present Data with Petukhov an- 
Popov, Brinich and Graham, a Calculated Least Squares Fit, 
and an Extrapolated Least Squares Fit. 



errors, in that the actual values are not given in their 
study and had to be taken from a plot of Stanton number 
versus Reynolds number. 

Disregarding the small differences between the results 
and correlations, the trend appears consistent. Heat trans- 
fer rates can be seen to increase both with increasing 
Reynolds number and curvature of the channel. The presence 
of Taylor-Gortler vortices in the laminar and transition 
regime seems to be the major contributing factor; whereas 
increased turbulence, possibly caused by the formation of 
these vortices, seems to dominate at higher Reynolds numbers. 
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VI . RECOMMENDATIONS 



The amount of experimental work remaining in the study 
of heat transfer in curved channels is significant. The 
present experimental apparatus has been used over a range 
of Reynolds numbers throughout the laminar, transition and 
into the turbulent regime. Additional experiments probing 
further into the turbulent regime should be conducted to 
better understand the heat transfer processes occurring 
there. This would require a more powerful compressor to 
generate the higher Reynolds numbers and a larger power 
supply to sufficiently drive the heated surfaces. Finally, 
increasing the number of thermocouples at least two-fold 
and the introduction of pressure transducers which could 
be connected directly to the data acquisition system would 
both prove to increase the accuracy and reliability of the 
experimental apparatus. 

To gain further insight into the effect of curvature 
on the rate of heat transfer, additional boundary conditions 
for the test should also be considered. An apparatus possess- 
ing a heated wall on each side of the test section might be 
one such consideration. 
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APPENDIX A: EXPERIMENTAL UNCERTAINTY 



The uncertainties for the major variable in the experi- 
ments, were calculated in accordance with the method des- 
cribed by S. Kline and F. McClintoch [Ref. 29] . The estimates 
of the uncertainty in the measured quantities were made con- 
servatively. As a result, there is considerable confidence 
in the uncertainties as calculated. 

The following equations were used to calculate the 
uncertainties : 




dp . 2 , . 2 

+ (1/4 (( - p a - - r ) + (^) )] 



1/2 



air 



2 



1/2 




) ) 




(4) dNu = Nu ( (^) + ( 

h 



dD 



D 



hd 



hd 
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The values obtained for the uncertainties for the values 



in the curved 


section test run at a Reynolds number of 



13000 are: 



Quantity 


Uncertainty 


A 


. 0348 


a pl 


. 0019 


A 

c 


. 0438 


c 

p 


. 0020 


D hd 


.0236 


h 


. 0058 


K 


. 0008 


K . 
air 


.0038 


Nu 


.0127 


m 


.0050 


Q air 


.0058 


Re hd 


. 0289 


T blk 


. 0004 


T. 

in 


. 0004 


T . 
out 


. 0004 


T 

WO 


.0002 


T -T . 

out in 


.0027 


Y 


. 0024 


AP 


.0500 


AT 


. 0005 


^air 


.0054 


p air 


.0001 
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The major source of uncertainty in the average Nusselt number 
is the uncertainty in the pressure readings from the mano- 
meter, which affects the mass flow rate calculation. In 
order to obtain a more accurate pressure measurement, it 
has been recommended that pressure transducers be used 
instead of vertical manometers in future studies. 



68 



REPRODUCED AT GOVERNMENT EXPENSE 



« 



APPENDIX B: SAMPLE CALCULATIONS 

Figure 23 below, shows the major heat transfer components 
for each of the test sections. The sample calculations that 
follow, demonstrate the methods used by the computer to calcu- 
late these components as well as the Reynold's number, average 
heat transfer coefficient, and average Nusselt number for each 
set of data. The sample calculations are for the curved sec- 
tion, but those for the straight section are similar. Reynolds 
number, Nusselt number, and Dean number are based on hydraulic 
diameter . 




Figure 23. Energy Balance in Straight Section. 
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A. SAMPLE CALCULATION DATA 



V H 


= 


55.28 V 


V R 


= 


4.53 V 


r pr 


= 


2.0262 Q 


P atm 


= 


30.11 in. Hg 


AP 


= 


12.00 in. H 2 0 


P 1 


= 


8.50 in. H^O 


T. 

m 


= 


23.01 C 


T . 
out 


= 


26.78 C 


T 

wo 


= 


44.74 C 


T . 

Wl 


= 


24.64 C 


T. n 
msl 


= 


39.67 C 


T . 

ms2 


= 


33.38 C 


T . . 
onf 


= 


26.57 C 


K . 
air 


= 


26.22 x 10' 3 W/mK 


C pair 


= 


1006 J/Kg K 


^air 


= 


18.59 x 10 ^ Kg/ms 


K. 

ms 


= 


3.6495 x 10” 2 W/mK 


AX . 
ms 


= 


0.01270 m 


£ 

wo 


= 


0.70 


£ • 
Wl 


= 


0.40 


a 


= 


5.669 x 10 -8 W/m 2 K 


S 


= 


. 5325 


Y 


= 


1.402 


gc 


= 


1 Kg m/NS 2 


R 


= 


286.8 Nm/Kg K 
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1.0 



wo-wi 



PL 



A 



A . 
pipe 

D 

c 



D . 
pipe 



D 



orif 



P 



wet 



.0730 ra 2 
.0016 m 2 
.002027 m 2 
.00635 ra 
.0508 ra 
.027051 m 
.5207 ra 



B. TEMPERATURE CALCULATIONS 

1. Bulk Temperature (T, , ) 

t blk 



T blk 



T. 



in 



+ T 



out 



2 



23.01 + 26.78 
2 



24.90 C 



2. Mean Temperature Difference (AT) 



AT = T -T, ... = 44.74 - 24.90 = 19.85 C 

wo blk 



3. Temperature Difference in the Insulation (AT. ) 
r ins 



AT. = T. . - T. 0 = 39.67 - 33.38 
ins insi ins2 



6.29 C 



C. POWER CALCULATIONS 

1. Power Supplied (Q ) 
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123.59 W 



Q 



P 



V H V R _ (55.28) (4.53) 

R pR (2.0262) 



2. Heat Lost Through the Outer Plate (Q 1q ) 



AT. 

ins 



-lo 



(AX. ) / (K. ) (A ) 

ins ins PL 



6.29 

(0.0127)/ (3.6495x10-^) (.0730) 



= 1.32 W 



3. Heat Radiated (Q ) 

a. Radiation Resistance (R R ) 



1 -e 



wo 



R n — ~ 

R e 

PL wo 



A F 
PL wo-wi 



1 -e . 
wi 

A PL £ wi 



i [_i_ + . i, . 



A„ T ' e e . 

PL wo wi 



[A + A-l]=40.12 m 2 



(0.0730) l 0.7 0.4 



b. Heat Radiated (Q r ) 



Qr = 



C (T 2 - T . 4 ) 
wo wi 



T [ = ] K 



= (5.669xl0~ 8 ) [ ( 44 . 74+273 ) 4 - ( 24 . 64+273 ) 4 ] 

40.12 



= 3.31 W 

D. MASS FLOW RATE CALCULATIONS 
1. Pressure Conversions 



a. P . =30.11 in. Hg x 3374.1 = 101594.15 N/m' 

atm ^ 
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